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CLUSTER EMISSION UNDER FEMTOSECOND LASER ABLATION OF SILICON 
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Rich populations of clusters have been observed after femtosecond laser ablation of bulk silicon 
in vacuum. Size and velocity distributions of the clusters as well Eis their charge states have been 
analyzed by reflectron time-of-flight mass spectrometry. An efficient emission of both neutral silicon 
clusters Sin (up to n = 6) and their cations Si^ (up to n = 10) has been observed. The clusters 
are formed even at very low laser fluences, below the ablation threshold, and their relative yield 
increases with fluence. We show the dependencies of the cluster yield as well as the expansion 
dynamics on both laser wavelength and laser fluence. The mechanisms of the cluster formation are 
discussed. 
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The recent development of femtosecond lasers expands 
the possible applications of laser ablation. However, the 
fundamental mechanisms of light-material interactions 
are still poorly understood. A study of the plume com- 
position and its expansion helps in understanding of the 
fundamental processes involved into the interaction and 
in the development of new laser applications. Recently, 
surface modifications under femtosecond laser irradiation 
of bulk silicon have been analyzed for laser pulse dura- 
tions from 5 to 400 fs 0|. The pump-and-probe studies 
show the moment of the plasma initiation and allow one 
to follow the plume expansion at the very early stages 
close to the target • The first observation of the silicon 
plume expansion at larger time scale after the femtosec- 
ond ablation by time-of-flight methods was reported in 
0, where the authors also mentioned the observation 
of small nanoclusters. However, very few information is 
available on the formation of silicon nanoclusters after 
the ablation by short laser pulses. The dimer desorption 
under resonant nanosecond excitation has been recently 
reported in the case of (2x1) reconstructed (100) silicon 
surface Q. The use of high energy photons of 6.4 eV 
(ArF laser) induces very strong desorption and ablation 
of atomic species and small clusters from both the bulk 
material and the nanostructured Si surfaces 0- How- 
ever, there is no systematical information neither on com- 
position nor on the dynamical properties of the plumes 
induced by femtosecond irradiation. In this work, we 
present evidence for the formation of clusters with sizes 
up to 10 atoms under ultrashort laser ablation of bulk 
silicon and analyze the dynamics of their expansion. 
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I. EXPERIMENTAL 

The experiments were performed with Si [100] surface 
under ultrahigh vacuum conditions 10^^° mbar). The 
monocrystalline Si target was irradiated at an angle of in- 
cidence of 45 ° using a Tiisapphire laser system (Mai-Tai 
coupled with TSA amplifier, Spectra Physics, 80 fs pulse 
duration, 10 Hz repetition rate, up to 30 mJ energy per 
pulse) operating at wavelengths of 800, 400, and 266 nm. 
A part of the laser beam was selected by an aperture to 
provide a nearly uniform intensity distribution over the 
irradiated spot. The target was rotated/translated dur- 
ing measurements to avoid considerable cratering. For 
each wavelength, series of craters at several laser ener- 
gies were produced on the stationary sample. By mea- 
suring the crater areas, the absolute calibration of laser 
fluence was performed. The laser fluence on the target 
was varied in the ranges 80—800, 20—200, 5—50 mJ/cm^ 
at wavelengths of 800, 400, and 266 nm, respectively. 

The abundance distributions of neutral and cationic 
particles in the laser-induced plume were analyzed by a 
reflectron time-of-flight mass spectrometer (TOF MS). 
When ions were studied, the plume was allowed to ex- 
pand under field- free conditions. While the neutral par- 
ticles of the plume were analyzed, the ionized species were 
rejected using a simple plasma suppressor j^. The latter 
consisted of a pair of deflection plates placed along the 
plume axis in front of the ion source of the MS where 
neutral particles were ionized by impact of 110 eV elec- 
trons. At a distance of 11 cm from the target, the ions, 
either ablated or post-ionized, were sampled parallel to 
the plume expansion axis by a 500 V repeller pulse af- 
ter a delay time td, following the laser pulse. All mass 
spectra were averaged over 300 laser shots. 



II. RESULTS AND DISCUSSION 

Efficient emission of both neutral silicon clusters Si„ 
(up to n = 6) and their cations Si^ (up to n ~ 10) have 
been observed for all investigated wavelengths. Si atoms 
and Si"*" ions were the most abundant particles at all laser 
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FIG. 1: Mass spectra of cationic (a) and neutral (b) species. 
The irradiation conditions are indicated in the figure. 



FIG. 2: Typical time-of-flight (TOF) distributions of Si+ ions 
for different fluences. Laser wavelength was 800 nm. Fth = 
100 mJ/cm^ is the threshold fluence corresponding to the ion 
appearance in the plume. 



fluences. Though, at certain irradiation regimes, the clus- 
ter fraction exceeded 10% of that of Si atoms and ions 
in the ablation plume. Fig. ^ shows typical mass spec- 
tra of cationic and neutral plume species, respectively. 
The spectra were obtained under the irradiation condi- 
tions corresponding to the maximum relative yield of the 
clusters with respect to monatomic particles. The abun- 
dance distributions are smooth under these conditions 
with peak intensities monotonously decreasing with clus- 
ter size. As it will be discussed below, the velocity dis- 
tributions for ionized plume particles are strongly depen- 
dent on particle mass. The ion mass spectra are , there- 
fore, different for particles expanding in the faster and the 
slower parts of the plume. The mass spectrum shown in 
Fig. ^ was obtained at the time delay corresponding to 
the maximum yield of Si^ dimer. 

A variation in the time delay allowed the analysis of 
particle velocity distributions and the characterization of 
the temporal evolution of the plume composition. Figure 
121 shows typical TOF distributions of Si"*" ions for differ- 
ent fluences. We have found that the distributions are 
weakly affected by both laser wavelength and fluence. In 
all cases, they reach the maximum at ~9 fis time delay 
that corresponds to the ion velocity of around 12 km/s, 



or to a kinetic energy of '^20 eV. An increase in laser flu- 
ence results only in a broadening of the distributions with 
a negligible shift towards higher velocities. At the same 
time, the total number of Si"*" ions increases strongly with 
fluence. The ion yield at a given fluence can be evalu- 
ated by integrating the TOF distributions in time. To 
obtain the values proportional to the total number of 
the plume particles at various fluences, this procedure 
should be corrected for angular distributions of expand- 
ing particles as a function of their velocity (usually fast 
plume species are more forward directed than slow ones 
The angular distributions were not measured in 
this work, but, since the ion velocity distributions are 
similar for different conditions (Fig|2), such a correction 
is not necessary for Si"*" ions. 

Figure |31 shows the yields of plume ions as a function 
of laser fluence for three laser wavelengths. These depen- 
dencies show a steep gradient that is nearly identical for 
800 and 400-nm wavelengths. For 266-nm, it is slightly 
weaker. We note that the total number of detected par- 
ticles varies by around 6 orders of magnitude in the stud- 
ied fluence ranges. The threshold fluences, Fth, for ion 
appearance in the plume, as deduced from MS measure- 
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FIG. 3: The ion yield as a function of laser fluence for 800, 
400 and 266 nm wavelengths of irradiation. 



mcnts arc around 100, 40, and 10 mJ/cm^ for 800, 400, 
and 266 nm, respectively. It should be noted that our ion 
detection system has a high sensitivity to detect a single 
plume ion. As a result, the obtained values correspond 
to the real thresholds for ion appearance but not to the 
ion detection threshold. 

Cluster ions Si^ were observed throughout the laser 
flucncc range studied, even at very low fluences close 
to the threshold for Si ions (at least for dimcrs). The 
clusters are found to be considerably slower than atomic 
ions. At near-threshold fluences, the TOF distributions 
for SiJ dimers are maximized at td =18 fj,s for all wave- 
lengths. The most probable velocity of the dimers is, 
therefore, approximately two times lower than that of 
Si^. This means that these plume particles have nearly 
the same momentum under these conditions. However, 
the Sig" trimers are faster than it would be expected for 
constant momentum (the TOF distributions arc maxi- 
mized at around 21 /is). 

By integrating the cluster TOF distributions in time, 
the total number of clusters in the plume as a function 
of fluence was also evaluated. For low fluences, the to- 
tal number of the detected dimers increases with flucncc 
and reaches maximum at F/Fth = 2 for 800 nm and 
at F/Fth — 4 for 400 nm. The maximum values of the 
population ratios SiJ/Si"'" are around 0.18 and 0.22 for 
800 and 400 nm respectively. Further increase in fluence 
results in an abrupt fall of the relative concentration of 
the dimers. This decrease of the SiJ abundance in the 
plume is due to both overall decrease of the relative clus- 



ter yield at high fluences and to a conversion of SiJ into 
larger cluster. The performed TOF integration shows the 
fast increase of the 814" and Si^ cluster abundance with 
fluence at F > 400 mJ/cm^. 

In contrast to Si"*" ions, the TOF distributions of clus- 
ters strongly depend on the laser fluence. Figure^jshows 
the evolution of the TOF distribution for the SiJ dimer 
at different fluences. These results were obtained for ex- 
citation at 800 nm. At very low fluences, below 200 
mJ/cm^ (= 2Fth), the distributions are rather narrow, 
single-peaked, and maximized at ~18 /xs. For laser flu- 
ences above 2Fth , the behavior of Si2^ changes dramat- 
ically. In the fluence range between 2Fth and 4:Ftin the 
TOF distributions are still single-peaked but their max- 
ima are significally shifted towards higher time delays 
(lower velocities). At F =4:Fth, the maximum is at '^25 
/iS (the corresponding most probable velocity is ~ 4.4 
km/s). The fast cluster ions are still present in this ab- 
lation regime, as can be seen in Fig. 01 but their distri- 
bution is masked by slower ions. With further increase 
in laser fluence the second slower population of Si^ ions 
with the most probable velocity of around 2.7 km/s ap- 
pears in the plume. At fluences above ^5.5Fth, this slow 
population becomes dominant. The first faster peak in 
the distribution (corresponding to the cluster ions de- 
tected in the (2— 4)_F't;i fluence range) is still present, well 
separated from the second peak, and is still maximized at 
~25 ^s. However, its intensity is much lower then at low 
fluences. The fast cluster ions (those observed at F < 
2Fth) are not detected anymore. 

The fluence dependency of the TOF distributions for 
larger Si^ clusters (n = 3—6) under under 800-nm ex- 
citation is found to be qualitatively similar to that for 
Si J. At a threshold fluence of around 450 mJ/cm^, the 
distributions are transformed from the single-peaked one 
to the double-peaked with the second (slow) cluster pop- 
ulation becoming rapidly dominant with further increase 
in flucncc. The most probable velocity of the second pop- 
ulation slightly decreases with cluster size from 2.5 km/s 
for Si J to 2 km/s for Si^. The abundance distribution for 
this slow cluster population is remarkably different from 
that for fast clusters (shown in Fig. ^). The Si4 and Sig 
cluster ions are more abundant than their odd-numbered 
neighbors. For larger time delays, the Sig cluster becomes 
dominant particle in the plume. 

Based on the obtained results, we can address the fun- 
damental questions on the mechanism of cluster forma- 
tion. Are the observed clusters formed primarily by their 
direct ejection from the target or by gas-phase aggre- 
gation of ablated atoms? Which process is responsible 
for such a transformation of the cluster velocity distribu- 
tion as shown in Fig. ^J" The set of the obtained TOF 
distributions and fluence dependencies provides a clear 
evidence for cluster formation mechanisms at different 
irradiation regimes. At very low laser fluence, from ion 
appearance threshold Fth up to ^2Fth, Si"*" and Si J ions 
are produced by an impulsive Coulomb explosion (CE) 
from a charged surface. We argue this by the following: 
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this field. The time of electrostatic interaction leading to 
the ion removal is thus many orders of magnitude shorter 
than the ion flight time to the mass spectrometer that re- 
sults in equal momenta for the ejected Si^ and SiJ ions as 
seen experimentally (Fig0}. We suggest that the electric 
field under these conditions is strong enough to repulse 
positive ions from the surface but not so strong to break 
the bonds of ejected clusters. Note that thermally des- 
orbed particles are characterized by equal kinetic energy 
rather than momentum 0. 

(b) An additional support for the CE is the high ve- 
locity of SiJ ions at low fluences. For fluences beyond 
^2Fth, the most probable velocity decreases that indi- 
cate the appearance of an additional ablation channel. 
For thermally desorbed particles one would expect the 
opposite behavior with fluence. 

Contrarily, the slower Si+ clusters (which overshadow 
the Coulomb explosion ions in the 2Fth—^Fth fluence 
range and form the fast cluster population at higher 
fluences) can be interpreted as "plasma ions", i.e., ions 
formed in the gas-phase ionized vapor plume. Also, we 
should take in consideration the direct thermal evapora- 
tion of small clusters from the surface (desorption of Sin 
clusters up to n = 6 was observed from thermally heated 
Si 0)- The velocities of these "plasma ions" decrease 
only slightly with cluster size and scale by a law inter- 
mediate between constant kinetic energy and constant 
velocity. 



FIG. 4: TOF distributions of Si2for irradiation with different 
fluences at laser wavelength of 800 nm. 



III. CONCLUSIONS 



(a) The observed momentum scaling for the particles 
and a very weak dependence of their expansion velocities 
on laser fluence and wavelength in the low fluence regimes 
indicates directly the CE mechanism 0. The repulsive 
electric field induced in the target by the fs-lascr pulse 
lasts for a short period of time (~1 ps). Both Si+ and 
SiJ spend roughly the same time in the action range of 



In conclusion, we present the first analysis of expan- 
sion dynamics of monoatomic species and small clusters 
after femtosecond laser ablation of clean (100) Si surface. 
While our results unambiguously show the domination 
of nonthermal mechanism of nanocluster emission at low 
laser fluence, more studies is needed to be certain about 
the origin of the ions under all the irradiation conditions. 
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